home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DPBFA.z / DPBFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPBBBBFFFFAAAA((((3333FFFF))))                                                            DDDDPPPPBBBBFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPBFA   - DPBFA factors a double precision symmetric positive definite
  10.      matrix stored in band form.
  11.  
  12.      DPBFA is usually called by DPBCO, but it can be called directly with a
  13.      saving in time if  RCOND  is not needed.
  14.  
  15.  
  16. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  17.       SUBROUTINE DPBFA(ABD,LDA,N,M,INFO)
  18.  
  19. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  20.      On Entry
  21.  
  22.      AAAABBBBDDDD DOUBLE PRECISION(LDA, N)
  23.         the matrix to be factored.  The columns of the upper
  24.         triangle are stored in the columns of ABD and the
  25.         diagonals of the upper triangle are stored in the
  26.         rows of ABD .  See the comments below for details.
  27.  
  28.      LLLLDDDDAAAA INTEGER
  29.         the leading dimension of the array  ABD .
  30.         LDA must be .GE. M + 1 .
  31.  
  32.      NNNN INTEGER
  33.         the order of the matrix  A .
  34.  
  35.      MMMM INTEGER
  36.         the number of diagonals above the main diagonal.
  37.         0 .LE. M .LT. N .  On Return
  38.  
  39.      AAAABBBBDDDD an upper triangular matrix  R , stored in band
  40.         form, so that  A = TRANS(R)*R .
  41.  
  42.      IIIINNNNFFFFOOOO INTEGER
  43.         = 0  for normal return.
  44.         = K  if the leading minor of order  K  is not
  45.         positive definite.  Band Storage
  46.         If  A  is a symmetric positive definite band matrix,
  47.         the following program segment will set up the input.
  48.         M = (band width above diagonal)
  49.         DO 20 J = 1, N
  50.         I1 = MAX0(1, J-M)
  51.         DO 10 I = I1, J
  52.         K = I-J+M+1
  53.         ABD(K,J) = A(I,J)
  54.         10    CONTINUE
  55.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  Cleve Moler,
  56.      University of New Mexico, Argonne National Lab.  Subroutines and
  57.      Functions BLAS DDOT Fortran MAX0,DSQRT
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.